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Abstract 

An analytical solution for symmetric Skyrmion was proposed for the SU(2) Skyrme model, which take the form of the hybrid form 
pf a kink-like solution and that given by the instanton method. The static properties of nucleons was then computed within the 
framework of collective quantization of the Skyrme model, with a good agreement with that given by the exact numeric solution. 
The comparisons with the previous results as well as the experimental values are also given. 
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1. Introduction 

The Skyrme model |Q]] is an effective field theory of mesons 
and baryons in which baryons arise as topological soliton so- 
lutions, known as Skyrmions. The model is based on the pre- 
QCD nonlinear <x model of the pion meson and was usually 
regarded to be consistent with the low-energy limit of large- 
N QCD|H. For this reason, among others, it has been exten- 
sively revisited in recent years [j, 13, S S §] (see, ifiol [TUl . 
for a review). Owing to the high nonlinearity, the solution to 
the Skyrme model was mainly studied through the numerical 
approach. It is worthwhile, however, to seek the analytic solu- 
tions IF2, 12l[l3l of Skyrmions due to its various applications 
in baryon phenomenology. One of noticeable analytic method 
for studying the Skyrmion solutions is the instanton approach 
proposed by Atiyah and Manton|8] which approximates criti- 
cal points of the Skyrme energy functional. 

In this Letter, we address the static solution of hedgehog 
Skyrmion in the S U (2) Skyrme model without pion mass term 
and propose an analytical solution for the hedgehog Skyrmion 
by writing it as the hybrid form of a kink-like solution and the 
analytic solution obtained by the instanton method [H]. Two 
lowest order of Pade approximations was used and the corre- 
sponding solutions for Skyrmion profile are given explicitly by 
using the downhill simplex method. The Skyrmion mass and 
static properties of nucleon as well as delta was computed and 
compared to the previous results. 

2. Analytic solution to Skyrme model 

The S U (2) Skyrme action without pion mass term is 
given by 
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d x 



Tr([L^L v f) 



(1) 



in which L M = U^d^U, U(x,t) e SU(2) is the nonlinear real- 
ization of the chiral field describing the <x field and n mesons 
with the unitary constrain U'U = 1, 2f„ the pion decay con- 
stant, and e a dimensionless constant characterizing nonlinear 
coupling. The Cauchy-Schwartz inequality for ([TJ implies lfl5ll 
E SK > 67T 2 (f K /e)\B\, where B = (1/24tt 2 ) j d*xs ijk Tr(LiLjL k ) 
is the topological charge, known as baryon number. Using the 
hedgehog ansatz, U(x) = cos(F) - i(x ■ <f) sin(F) (cr are the three 
Pauli matrices) with F = F(r) depending merely on the radial 
coordinate r, the static energy for ([1) becomes 



7 SK 



= 2tt( 



¥1 



dx 



x 2 F 2 x + 2sin z (F)(l +F l x ) + 



sin 4 (F) 



(2) 

with x - ef„r a dimensionless variable and F x = dF(x)/dx. 
The equation of motion of (O is 



sin 2 F\ 2 sin(2F). , 



sin 2 F \ 



= 0, (3) 



where the boundary condition F(Q) = n, F(oo) = will be 
imposed so that it corresponds to the physical vacuum for 
(Q.U = +1. The equation (01 is usually solved numerically 
due to its high nonlinearity yl |4j, |6[ |5[ . 

A kink-like analytic solution was given by |5tl 



Fi(x) = 4 arctan[exp(-x)], 



(4) 



with E SK = l.24035(6n 2 f„/e), while an alternative Skyrmion 
profile, proposed based on the instanton method, takes the 
formli 

A 

y2' 



F 2 (x) = n 



i-(i + -rr 1/2 



(5) 
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with corresponding energy 1.2432(6:7r 2 / 7r /e) for the numeric 
factor A = 2.109. The singularity at r = is gauge depen- 
dent and can be gauged away without affecting the value for the 
Skyrme field. 

To find the more accurate analytic solution, we first improve 
the solution (0 into 4 arctan[exp(-cx)] with c a numeric fac- 
tor and then take A to be a x-dependent function: A — » A(x). 
Hence, we propose a Skyrmion profile function in the hybrid 
form mixing (0 and (0 



F w (x) = 4warctan[exp(-cx)] + n{\ - w) 



i-<i + S- 



1/2 



(6) 

with w € [0, 1] being a positive weight factor. In principle, 
one can find the governing equation for the unknown A(x) by 
substituting (0 into and obtain a series solution of A(x) by 
solving the governing equation. Here, however, we choose the 
Pade approximation to parameterize A(x) 



A(x) = A 



1 + ax 2 + ■ ■ ■ 
1 + bx 2 + ■ ■ • ' 



(7) 



since it has as equal potential as series in approximating a con- 
tinuous function. Note that we have already written A(x) as 
function of x 2 instead of x since so is f^M in (0. The simplest 
nontrivial case of the above Pade approximation is the [2/2] ap- 
proximant 

1 + ax 2 



A(x) = /Iq- 



1 + bx 2 ' 

The minimization of the energy (0 with the trial function 
(0 with respect to the variational parameters (a, b, c, w, Aq) was 
carried out numerically for the [2/2] Pade approximant ([8]) us- 
ing the downhill simplex method (the Neilder-Mead algorithm). 
The result for the optimized parameters is given by 



a = 0.330218, b = 1.331975, c= 2.094056, 
w = 0.286566, ^ = 7.323877. 



(9) 



with E SK = 1.23152(67r 2 /^/e)- The solution ©, with A(x) 
given by (0 and the parameters (0, is referred as solution 
Hyb(2/2) for short in this paper and is plotted in Fig.l, com- 
pared to the solutions (0 and (0, and the numerical solution 
(Num.) to the equation (0. We also include the analytic solu- 
tions given by [14] and the solution in the form of the purely 
Pade approximant | for comparison. A quite well agreement 
of our solution with the numerical solution can be seen from 
this plot. We note that the inequality E s K / (6n 2 f K / e) > \B\ = 1 
is fulfilled for all of cited results of the energy(see Table I). 

To check how well the asymptotic behavior of (0 is we apply 
the asymptotic expansion analysis on the profile F(x). For small 
x <k 1 the solution F(x) to the equation (0 is given by 



F(x) =jt + Ax + Bx 3 +Cx 5 + ■■■ 

=tt - 2.007528x + 0.358987x 3 - 0.146499x 5 + 



(10) 



(see also lll2ll . where the variable x used is twice of x in this 
paper) while the analytic solution (0, when (0 and (0 is used, 



behaves like 

F w (x) = tt-2.028368x + 0.518855x 3 -0.641539x 5 + --- . (11) 

One can see that (Q~TJ| agrees well with (TTOb up to or. For large 
x — > oo the series solution for F(x) can be obtained by solving 
(0 with x replaced by 1 /y and using the series expansion for 
small y. After re-changing y to x, one finds 

, s 2.1596 f 0.222 116.0 0.113 2.71 

F(X) = — I 1 ; — + — =- + — + • 



(12) 

On the other hand, the solution Hyb(2/2), at large x, has the 
asymptotic form 



2.0348 f 0.91576 5.8524 9.6819 
F w (x) = =— 1 1 + 5 — + 



3.1677 51.943 



(13) 



which agrees globally with ( fT2T > except for a small bit differ- 
ences. The detailed differences between ( fT2b and ( TT3T > at large x 
can be due to the fact that the variationally-obtained solution (0 
approximates the Skyrmion profile globally and may produce 
small errors in local region, for instance, F W (5Q) = 8. 142 x 10~ 4 
while F(50) = 8.638 x 10~ 4 . 

The disagreement can be improved by employing the Pade 
approximant of higher order than (0, for example, the [4/4] 
approximant 

1 + ax 2 + diX A 
A{x) - A \ + bx 2 + b^ (U) 

The minimization of (0 using (TBI) , as done for the [2/2] ap- 
proximant, yields the numerically optimal parameters, 



a = 0.2598, b = 0.5446, c = 1 .9932, a 4 = 0.0538, 
b A = 0.1226, w = 0.1839, A = 3.9439, 



(15) 



The solution (0 with A(x) specified by ( TBI and ( TT3T > will be re- 
ferred as the Hyb(4/4) in this paper and is also plotted in Fig.l. 
The Fig. 2 shows the profiles of F(x) at large x for Hyb(2/2) as 
well as Hyb(4/4), and the numeric solution. The asymptotic ex- 
pansion of the solution Hyb(4/4) shows that for small x <K 1 the 
profile becomes 

FJx) =n- 2.0243x + 0.4654x 3 - 0.4270x 5 + • • • . 



while for x — » 

F w (x) 



2.220 



0.9150 9.5908 65.358 



x L 

264.78 



of 

821.88 



Here, a better value F w (50) = 8.877 x 10~ 4 is obtained for the 
latter asymptotic profile in contrast with the solution Hyb(2/2). 
The computed Skyrmion energies (0, measured in the unit of 
2f w fe, are listed in Table I, including the corresponding results 
obtained by the numeric solution and obtained in the relevant 
references as indicated. 
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In solving (01 numerically, we employ the nonlinear shoot 
algorithm for the boundary values at x — 0.001 and x = 1000 
based on the asymptotic formulas ( ITOb and ( fT2T > of the chiral 
angle F(x). 

3. The static properties of nucleons at low energy 

The static properties of nucleons can be extracted by semi- 
classically quantizing the spinning modes of Skyrme La- 
grangian using the collective variables yy. Here, we will use 
the solution Hyb(2/2) and Hyb(4/4) to compute the static prop- 
erties of nucleons and nucleon-isobar(A) within the framework 
of the bosonic quantization of Skyrme model. 

Following Adkin et al.j^], one can choose S t/(2)-variable 
A(t) as the collective variables, and substitute U — 
A(t)Uo(x)A(t) Jf into (|TJ. In the adiabatic limit, one has 



S SK = S S 



dtTr\ 



r dA <9A + 
~dt~dT 



(16) 



with Sq K the action for the static hedgehog configuration, Iq 

n/(.3e 3 f„), and 



Jr»oo 
x : 
o 



dxsin 2 F[l + F 2 X + sin 2 F/x 2 ]. 



(17) 



which is independent of f„ and e. The Hamiltonian associ- 
ated to (fTol l. when quantized via the quantization procedure 
in terms of collective coordinates, yields an eigenvalue (H) = 
M + J(J + l)/(27oA), with M - E SK being the soli ton energy 
of the Skyrmion. This yields the masses of the nucleon and 
A-isobar 

M /v=M+-^-,Ma = M+-^-. (18) 
8/oA 8/ ( )A 

By adjusting f n and e to fit the nucleon and delta masses in ( fT8l . 
one can fix the model parameters f„ and e using the calculated 
M and A through © and ( TPTT i. 

The isoscalor root mean square(r.m.s) radius and isoscalor 
magnetic r.m.s radius are given by 



efn{r 2 ))L 2 ={— x 2 sin 2 FF X 



1/2 



efir 2 ) 1 ' 2 = 



C x A sin 2 FF x dx 1 
C° x 2 sin 2 FF x dx j 



1/2 



(19) 



respectively. Combining with the masses of nucleon and the 
delta, one can evaluate the magnetic moments for proton and 
neutron via the following formula 



flp,n 



<r 2 >/= 



M N (M A - M N ) : 



2{M A -M N Y 
(20) 

where plus and minus correspond to proton and neutron, respec- 
tively. The calculated results for these quantities using two so- 
lution schemes (Hyb.(2/2) and Hyb.(4/4)) are shown explicitly 
in Table II, compared to the experimental values as well as that 
computed by the numeric solution for F(x). The corresponding 
results from other predictions are also shown in this table. Here 



in Table II, we use the experimental values = 938.9MeV, 
Ma = 1232MeV for fixing e and f n through ( fTSl . in contrast 
with me input Ma, = 938MeV,M A = 1232MeV usedby Ref.Q 
and Ref.lfll. 

To check the solution further, we also list, in the Table II, the 
axial coupling constant and the 7rAW-coupling, which are given 
by 

gA = —^G,g„ NN = —g A (21) 



respectively. Here, the numeric factor G is 



Jo 



4 I dxxr 
Jo 



F x + 



sin 2F sin 2F 



2 sin 2 F sin 2 F 
r — F x + — r— sm2F 



(22) 



4. Concluding remarks 

We show that the hybrid form of a kink-like solution and 
that given by the instanton method are suited to approximate 
the exact solution for the hedgehog Skyrmion, when combining 
with Pade approximation. The resulted analytic solution (0 
has two remarkable features: (1) it is simple in the sense that 
it is globally given in whole region; (2) it well approaches the 
asymptotic behavior of the exact solution. We note that the 
further generalization of ©, made by approximating c in © via 
Pade approximation, does not exhibit remarkable improvement, 
particularly in the asymptotic behavior of the chiral angle F (x) 
at infinity. We expect that our solution can be useful in the 
dynamics study of the Skyrmion evolution and interactions. 
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Table I 



Work 


m 


[13] 


16] 


[14] 


[12] 


Hyb(2/2) 


Hyb(4/4) 


Num. 


(2/,/e) 


36.5 


36.47 


36.484 


36.474 


36.47 


36.4638 


36.4638 


36.4616 


M 


1.233 


1.2317 


1.2322 


1.23186 


1.2317 


1.23152 


1.23146 


1.23145 



Table II 



Quantities 


Ref.[3] 


Ref.[14] 


Hyb.(2/2) 


Hyb.(4/4) 


Num. 


Expt. 


M/(2fJe) 


36.5 


36.5 


36.4638 


36.4638 


36.4616 




2f n (MeV) 


129 


130 


128.730 


129.453 


129.260 


186 


e 


5.45 


5.48 


5.4229 


5.4527 


5.4446 




A 


50.9 


52.2 


50.1467 


51.2830 


50.9782 




(r 2 )\'i{fm) 


0.59 


0.586 


0.5985 


0.5920 


0.5938 


0.72 




0.92 


0.920 


0.9258 


0.9208 


0.9222 


0.81 




1.87 




1.8825 


1.8764 


1.8781 


2.79 




-1.31 




-1.3209 


-1.3269 


-1.3253 


-1.91 




1.43 




1.4252 


1.4141 


1.4171 


1.46 


8a 


0.61 




0.6332 


0.5992 


0.6081 


1.23 


gnNN 


8.9 




9.2364 


8.6921 


8.8343 


13.5 
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1. Introduction 

The Skyrme model |Q]] is an effective field theory of mesons 
and baryons in which baryons arise as topological soliton so- 
lutions, known as Skyrmions. The model is based on the pre- 
QCD nonlinear <x model of the pion meson and was usually 
regarded to be consistent with the low-energy limit of large- 
N QCD|H. For this reason, among others, it has been exten- 
sively revisited in recent years [j, 13, S S §] (see, ifiol [TUl . 
for a review). Owing to the high nonlinearity, the solution to 
the Skyrme model was mainly studied through the numerical 
approach. It is worthwhile, however, to seek the analytic solu- 
tions IF2, 12l[l3l of Skyrmions due to its various applications 
in baryon phenomenology. One of noticeable analytic method 
for studying the Skyrmion solutions is the instanton approach 
proposed by Atiyah and Manton|8] which approximates criti- 
cal points of the Skyrme energy functional. 

In this Letter, we address the static solution of hedgehog 
Skyrmion in the S U (2) Skyrme model without pion mass term 
and propose an analytical solution for the hedgehog Skyrmion 
by writing it as the hybrid form of a kink-like solution and the 
analytic solution obtained by the instanton method [H]. Two 
lowest order of Pade approximations was used and the corre- 
sponding solutions for Skyrmion profile are given explicitly by 
using the downhill simplex method. The Skyrmion mass and 
static properties of nucleon as well as delta was computed and 
compared to the previous results. 

2. Analytic solution to Skyrme model 

The S U (2) Skyrme action without pion mass term is 
given by 
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d x 



Tr([L^L v f) 



(1) 



in which L M = U^d^U, U(x,t) e SU(2) is the nonlinear real- 
ization of the chiral field describing the <x field and n mesons 
with the unitary constrain U'U = 1, 2f„ the pion decay con- 
stant, and e a dimensionless constant characterizing nonlinear 
coupling. The Cauchy-Schwartz inequality for ([TJ implies lfl5ll 
E SK > 67T 2 (f K /e)\B\, where B = (1/24tt 2 ) j d*xs ijk Tr(LiLjL k ) 
is the topological charge, known as baryon number. Using the 
hedgehog ansatz, U(x) = cos(F) - i(x ■ <f) sin(F) (cr are the three 
Pauli matrices) with F = F(r) depending merely on the radial 
coordinate r, the static energy for ([1) becomes 



7 SK 



= 2tt( 



¥1 



dx 



x 2 F 2 x + 2sin z (F)(l +F l x ) + 



sin 4 (F) 



(2) 

with x - ef„r a dimensionless variable and F x = dF(x)/dx. 
The equation of motion of (O is 



sin 2 F\ 2 sin(2F). , 



sin 2 F \ 



= 0, (3) 



where the boundary condition F(Q) = n, F(oo) = will be 
imposed so that it corresponds to the physical vacuum for 
(Q.U = +1. The equation (01 is usually solved numerically 
due to its high nonlinearity yl |4j, |6[ |5[ . 

A kink-like analytic solution was given by |5tl 



Fi(x) = 4 arctan[exp(-x)], 



(4) 



with E SK = l.24035(6n 2 f„/e), while an alternative Skyrmion 
profile, proposed based on the instanton method, takes the 
formli 

A 

y2' 



F 2 (x) = n 



i-(i + -rr 1/2 



(5) 



Preprint submitted to Phys lett B 



December 28, 2009 



with corresponding energy 1.2432(6:7r 2 / 7r /e) for the numeric 
factor A = 2.109. The singularity at r = is gauge depen- 
dent and can be gauged away without affecting the value for the 
Skyrme field. 

To find the more accurate analytic solution, we first improve 
the solution (0 into 4 arctan[exp(-cx)] with c a numeric fac- 
tor and then take A to be a x-dependent function: A — » A(x). 
Hence, we propose a Skyrmion profile function in the hybrid 
form mixing (0 and (0 



F w (x) = 4warctan[exp(-cx)] + n{\ - w) 



i-<i + S- 



1/2 



(6) 

with w € [0, 1] being a positive weight factor. In principle, 
one can find the governing equation for the unknown A(x) by 
substituting (0 into and obtain a series solution of A(x) by 
solving the governing equation. Here, however, we choose the 
Pade approximation to parameterize A(x) 



A(x) = A 



1 + ax 2 + ■ ■ ■ 
1 + bx 2 + ■ ■ • ' 



(7) 



since it has as equal potential as series in approximating a con- 
tinuous function. Note that we have already written A(x) as 
function of x 2 instead of x since so is f^M in (0. The simplest 
nontrivial case of the above Pade approximation is the [2/2] ap- 
proximant 

1 + ax 2 



A(x) = /Iq- 



1 + bx 2 ' 

The minimization of the energy (0 with the trial function 
(0 with respect to the variational parameters (a, b, c, w, Aq) was 
carried out numerically for the [2/2] Pade approximant ([8]) us- 
ing the downhill simplex method (the Neilder-Mead algorithm). 
The result for the optimized parameters is given by 



a = 0.330218, b = 1.331975, c= 2.094056, 
w = 0.286566, ^ = 7.323877. 



(9) 



with E SK = 1.23152(67r 2 /^/e)- The solution ©, with A(x) 
given by (0 and the parameters (0, is referred as solution 
Hyb(2/2) for short in this paper and is plotted in Fig.l, com- 
pared to the solutions (0 and (0, and the numerical solution 
(Num.) to the equation (0. We also include the analytic solu- 
tions given by [14] and the solution in the form of the purely 
Pade approximant | for comparison. A quite well agreement 
of our solution with the numerical solution can be seen from 
this plot. We note that the inequality E s K / (6n 2 f K / e) > \B\ = 1 
is fulfilled for all of cited results of the energy(see Table I). 

To check how well the asymptotic behavior of (0 is we apply 
the asymptotic expansion analysis on the profile F(x). For small 
x <k 1 the solution F(x) to the equation (0 is given by 



F(x) =jt + Ax + Bx 3 +Cx 5 + ■■■ 

=tt - 2.007528x + 0.358987x 3 - 0.146499x 5 + 



(10) 



(see also lll2ll . where the variable x used is twice of x in this 
paper) while the analytic solution (0, when (0 and (0 is used, 



behaves like 

F w (x) = tt-2.028368x + 0.518855x 3 -0.641539x 5 + --- . (11) 

One can see that (Q~TJ| agrees well with (TTOb up to or. For large 
x — > oo the series solution for F(x) can be obtained by solving 
(0 with x replaced by 1 /y and using the series expansion for 
small y. After re-changing y to x, one finds 

, s 2.1596 f 0.222 116.0 0.113 2.71 

F(X) = — I 1 ; — + — =- + — + • 



(12) 

On the other hand, the solution Hyb(2/2), at large x, has the 
asymptotic form 



2.0348 f 0.91576 5.8524 9.6819 
F w (x) = =— 1 1 + 5 — + 



3.1677 51.943 



(13) 



which agrees globally with ( fT2T > except for a small bit differ- 
ences. The detailed differences between ( fT2b and ( TT3T > at large x 
can be due to the fact that the variationally-obtained solution (0 
approximates the Skyrmion profile globally and may produce 
small errors in local region, for instance, F W (5Q) = 8. 142 x 10~ 4 
while F(50) = 8.638 x 10~ 4 . 

The disagreement can be improved by employing the Pade 
approximant of higher order than (0, for example, the [4/4] 
approximant 

1 + ax 2 + diX A 
A{x) - A \ + bx 2 + b^ (U) 

The minimization of (0 using (TBI) , as done for the [2/2] ap- 
proximant, yields the numerically optimal parameters, 



a = 0.2598, b = 0.5446, c = 1 .9932, a 4 = 0.0538, 
b A = 0.1226, w = 0.1839, A = 3.9439, 



(15) 



The solution (0 with A(x) specified by ( TBI and ( TT3T > will be re- 
ferred as the Hyb(4/4) in this paper and is also plotted in Fig.l. 
The Fig. 2 shows the profiles of F(x) at large x for Hyb(2/2) as 
well as Hyb(4/4), and the numeric solution. The asymptotic ex- 
pansion of the solution Hyb(4/4) shows that for small x <K 1 the 
profile becomes 

FJx) =n- 2.0243x + 0.4654x 3 - 0.4270x 5 + • • • . 



while for x — » 

F w (x) 



2.220 



0.9150 9.5908 65.358 



x L 

264.78 



of 

821.88 



Here, a better value F w (50) = 8.877 x 10~ 4 is obtained for the 
latter asymptotic profile in contrast with the solution Hyb(2/2). 
The computed Skyrmion energies (0, measured in the unit of 
2f w fe, are listed in Table I, including the corresponding results 
obtained by the numeric solution and obtained in the relevant 
references as indicated. 
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In solving (01 numerically, we employ the nonlinear shoot 
algorithm for the boundary values at x — 0.001 and x = 1000 
based on the asymptotic formulas ( ITOb and ( fT2T > of the chiral 
angle F(x). 

3. The static properties of nucleons at low energy 

The static properties of nucleons can be extracted by semi- 
classically quantizing the spinning modes of Skyrme La- 
grangian using the collective variables yy. Here, we will use 
the solution Hyb(2/2) and Hyb(4/4) to compute the static prop- 
erties of nucleons and nucleon-isobar(A) within the framework 
of the bosonic quantization of Skyrme model. 

Following Adkin et al.j^], one can choose S t/(2)-variable 
A(t) as the collective variables, and substitute U — 
A(t)Uo(x)A(t) Jf into (|TJ. In the adiabatic limit, one has 



S SK = S S 



dtTr\ 



r dA <9A + 
~dt~dT 



(16) 



with Sq K the action for the static hedgehog configuration, Iq 

n/(.3e 3 f„), and 



Jr»oo 
x : 
o 



dxsin 2 F[l + F 2 X + sin 2 F/x 2 ]. 



(17) 



which is independent of f„ and e. The Hamiltonian associ- 
ated to (fTol l. when quantized via the quantization procedure 
in terms of collective coordinates, yields an eigenvalue (H) = 
M + J(J + l)/(27oA), with M - E SK being the soli ton energy 
of the Skyrmion. This yields the masses of the nucleon and 
A-isobar 

M /v=M+-^-,Ma = M+-^-. (18) 
8/oA 8/ ( )A 

By adjusting f n and e to fit the nucleon and delta masses in ( fT8l . 
one can fix the model parameters f„ and e using the calculated 
M and A through © and ( TPTT i. 

The isoscalor root mean square(r.m.s) radius and isoscalor 
magnetic r.m.s radius are given by 



efn{r 2 ))L 2 ={— x 2 sin 2 FF X 



1/2 



efir 2 ) 1 ' 2 = 



C x A sin 2 FF x dx 1 
C° x 2 sin 2 FF x dx j 



1/2 



(19) 



respectively. Combining with the masses of nucleon and the 
delta, one can evaluate the magnetic moments for proton and 
neutron via the following formula 



flp,n 



<r 2 >/= 



M N (M A - M N ) : 



2{M A -M N Y 
(20) 

where plus and minus correspond to proton and neutron, respec- 
tively. The calculated results for these quantities using two so- 
lution schemes (Hyb.(2/2) and Hyb.(4/4)) are shown explicitly 
in Table II, compared to the experimental values as well as that 
computed by the numeric solution for F(x). The corresponding 
results from other predictions are also shown in this table. Here 



in Table II, we use the experimental values = 938.9MeV, 
Ma = 1232MeV for fixing e and f n through ( fTSl . in contrast 
with me input Ma, = 938MeV,M A = 1232MeV usedby Ref.Q 
and Ref.lfll. 

To check the solution further, we also list, in the Table II, the 
axial coupling constant and the 7rAW-coupling, which are given 
by 

gA = —^G,g„ NN = —g A (21) 



respectively. Here, the numeric factor G is 



Jo 



4 I dxxr 
Jo 



F x + 



sin 2F sin 2F 



2 sin 2 F sin 2 F 
r — F x + — r— sm2F 



(22) 



4. Concluding remarks 

We show that the hybrid form of a kink-like solution and 
that given by the instanton method are suited to approximate 
the exact solution for the hedgehog Skyrmion, when combining 
with Pade approximation. The resulted analytic solution (0 
has two remarkable features: (1) it is simple in the sense that 
it is globally given in whole region; (2) it well approaches the 
asymptotic behavior of the exact solution. We note that the 
further generalization of ©, made by approximating c in © via 
Pade approximation, does not exhibit remarkable improvement, 
particularly in the asymptotic behavior of the chiral angle F (x) 
at infinity. We expect that our solution can be useful in the 
dynamics study of the Skyrmion evolution and interactions. 
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Table I 



Work 


m 


[13] 


16] 


[14] 


[12] 


Hyb(2/2) 


Hyb(4/4) 


Num. 


(2/,/e) 


36.5 


36.47 


36.484 


36.474 


36.47 


36.4638 


36.4638 


36.4616 


M 


1.233 


1.2317 


1.2322 


1.23186 


1.2317 


1.23152 


1.23146 


1.23145 



Table II 



Quantities 


Ref.[3] 


Ref.[14] 


Hyb.(2/2) 


Hyb.(4/4) 


Num. 


Expt. 


M/(2fJe) 


36.5 


36.5 


36.4638 


36.4638 


36.4616 




2f n (MeV) 


129 


130 


128.730 


129.453 


129.260 


186 


e 


5.45 


5.48 


5.4229 


5.4527 


5.4446 




A 


50.9 


52.2 


50.1467 


51.2830 


50.9782 




(r 2 )\'i{fm) 


0.59 


0.586 


0.5985 


0.5920 


0.5938 


0.72 




0.92 


0.920 


0.9258 


0.9208 


0.9222 


0.81 




1.87 




1.8825 


1.8764 


1.8781 


2.79 




-1.31 




-1.3209 


-1.3269 


-1.3253 


-1.91 




1.43 




1.4252 


1.4141 


1.4171 


1.46 


8a 


0.61 




0.6332 


0.5992 


0.6081 


1.23 


gnNN 


8.9 




9.2364 


8.6921 


8.8343 


13.5 
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Abstract 

An analytical solution for symmetric Skyrmion was proposed for the SU(2) Skyrme model, which take the form of the hybrid form 
pf a kink-like solution and that given by the instanton method. The static properties of nucleons was then computed within the 
framework of collective quantization of the Skyrme model, with a good agreement with that given by the exact numeric solution. 
The comparisons with the previous results as well as the experimental values are also given. 
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PACS: 12.38.-t, 11.15.Tk, 12.38.Aw 
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1. Introduction 

The Skyrme model |Q]] is an effective field theory of mesons 
and baryons in which baryons arise as topological soliton so- 
lutions, known as Skyrmions. The model is based on the pre- 
QCD nonlinear <x model of the pion meson and was usually 
regarded to be consistent with the low-energy limit of large- 
N QCD|H. For this reason, among others, it has been exten- 
sively revisited in recent years [j, 13, S S §] (see, ifiol [TUl . 
for a review). Owing to the high nonlinearity, the solution to 
the Skyrme model was mainly studied through the numerical 
approach. It is worthwhile, however, to seek the analytic solu- 
tions IF2, 12l[l3l of Skyrmions due to its various applications 
in baryon phenomenology. One of noticeable analytic method 
for studying the Skyrmion solutions is the instanton approach 
proposed by Atiyah and Manton|8] which approximates criti- 
cal points of the Skyrme energy functional. 

In this Letter, we address the static solution of hedgehog 
Skyrmion in the S U (2) Skyrme model without pion mass term 
and propose an analytical solution for the hedgehog Skyrmion 
by writing it as the hybrid form of a kink-like solution and the 
analytic solution obtained by the instanton method [H]. Two 
lowest order of Pade approximations was used and the corre- 
sponding solutions for Skyrmion profile are given explicitly by 
using the downhill simplex method. The Skyrmion mass and 
static properties of nucleon as well as delta was computed and 
compared to the previous results. 

2. Analytic solution to Skyrme model 

The S U (2) Skyrme action without pion mass term is 
given by 
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in which L,, 



Tr(L II ) 2 + —Tr([L p ,L v ] 2 )] (1) 



U(x, t) e S U(2) is the nonlinear real- 



ization of the chiral field describing the cr field and n mesons 
with the unitary constrain U'U = 1, 2f„ the pion decay con- 
stant, and e a dimensionless constant characterizing nonlinear 
coupling. The Cauchy-Schwartz inequality for ([TJ, implies lfl5ll 
E SK > 67T Z (f K /e)\B\, where B = (1/24tt 2 ) j d 3 xj jk Tr(LiLjL k ) 
is the topological charge, known as baryon number. Using the 
hedgehog ansatz, U(x) = cos(F) - i(x ■ <f) sin(F) (cr are the three 
Pauli matrices) with F = F(r) depending merely on the radial 
coordinate r, the static energy for (Q]i becomes 



E SK = 2n( 



¥1 



dx 



->-> -> ? sin (F) 
x 2 F 2 x + 2 sin 2 (F)(l + F 1 ) + )-- 



(2) 

with x = ef„r a dimensionless variable and F x = dF(x)/dx. 
The equation of motion of (f2]i is 

/ sin 2 F\ 2 sin(2F) / , sin 2 F\ 

where the boundary condition F(0) — n, F (oo) = Q will be im- 
posed so that it corresponds to the physical vacuum for (|TJ: f/ = 
+ 1. The equation (0) is usually solved numerically due to its 
highly-nonlinearity, as done in most calculations yl|4j,|6t|5[]. 
A kink-like analytic solution was given by 



Fi(x) = 4arctan[exp(-x)], 



(4) 



with E SK = l.24035(6n 2 f„/e), while an alternative Skyrmion 
profile, proposed based on the instanton method, takes the 
formli 

A 



F 2 (x) = 7T 



1 - (1 + — ) 

x z 



■1/2 



(5) 

December 28, 2009 



with corresponding energy 1.2432(6:7r 2 / 7r /e) for the numeric 
factor A = 2.109. The singularity at r = is gauge depen- 
dent and can be gauged away without affecting the value for the 
Skyrme field. 

To find the more accurate analytic solution, we first improve 
the solution (0 into 4 arctan[exp(-cx)] with c a numeric fac- 
tor and then take A to be a x-dependent function: A — » A(x). 
Hence, we propose a Skyrmion profile function in the hybrid 
form mixing (0 and (0 



F w (x) = 4warctan[exp(-cx)] + n(l - w) 



(6) 

with w e [0, 1] being a positive weight factor. In principle, 
one can find the governing equation for the unknown A(x) by 
substituting (0 into and obtain a series solution of A(x) by 
solving the governing equation. Here, however, we choose the 
Pade approximation to parameterize A(x) 



A(x) = Aq 



1 + ax 2 H 

1 + bx 1 + • • • ' 



(7) 



since it has as equal potential as series in approximating a con- 
tinuous function. Note that we have already written A(x) as 
function of x 2 instead of x since so is F2(x) in (0. The simplest 
nontrivial case of the above Pade approximation is the [2/2] ap- 
proximant 

1 + ax 2 



A(x) = /Iq- 



1 + bx 2 ' 

The minimization of the energy (0 with the trial function 
(0 with respect to the variational parameters (a, b, c, w, Aq) was 
carried out for the [2/2] Pade approximant (0 using the down- 
hill simplex method (the Neilder-Mead algorithm). The result 
for the numerically optimized parameters is given by 



a = 0.330218, £ = 1.331975, c= 2.094056, 
w = 0.286566, Ao = 7.323877. 



(9) 



with E SK = 1.23152(67r 2 /^/e). The solution (0, with A(x) 
given by (0 and the parameters (0, is referred as solution 
Hyb(2/2) for short in this paper and is plotted in Fig.l, com- 
pared to the solutions (0 and @, and the numerical solution 
(Num.) to the equation 0. We also include the analytic solu- 
tions given by 1 1411 an d the solution in the form of the purely 
Pade approximant 111 311 for comparison. A quite well agreement 
of our solution with the numerical solution can be seen from 
this plot. We note that the inequality E s K / \6n 2 f„ / e) > \B\ = 1 
is fulfilled for all of cited results of the energy(see Table I). 

To check how well the asymptotic behavior of (0 is we apply 
the asymptotic expansion analysis on the profile F(x). For small 
x <3C 1 the solution F(x) to the equation (0 can be given by 
solving the equation via series expansion, which is 

F(x) =n + Ax + fix 3 + Cx 5 + ■■■ 

=n - 2.007528x + 0.358987x 3 - 0.146499x 5 + • • • , 

(10) 

(see also [12], where the variable x used is twice of x in this 
paper) while the analytic solution (0, when (0 is used, behaves 



like 

F w (x) = tt-2.028368x + 0.518855x 3 -0.641539x 5 + --- . (11) 

One can see that (fTTT i agrees well with (TTOb up to x . For large 
x — > oo the series solution for F(x) can be obtained by solving 
with x replaced by 1 /y and using the series expansion for 
small y. After re-changing y to x, one finds 

, s 2.1596 f 0.222 116.0 0.113 2.71 

F(X) = — I 1 ; — + — =- + — + • 



(12) 

On the other hand, the solution Hyb(2/2), at large x, has the 
asymptotic form 



2.0348 f 0.91576 5.8524 9.6819 
F w {x) = t— { 1 + — + 



3.1677 51.943 



(13) 



which agrees globally with $1% except for a small bit differ- 
ences. The detailed differences between ( fT2b and ( TT3T > at large x 
can be due to the fact that the variationally-obtained solution (0 
approximates the Skyrmion profile globally and may produce 
small errors in local region, for instance, F W (5Q) = 8. 142 x 10~ 4 
while F(50) = 8.638 x 10~ 4 . 

The disagreement can be improved by employing the Pade 
approximant of higher order than (0, for example, the [4/4] 
approximant 

1 + ax 2 + diX A 
A{x) - A \ + bx 2 + b^ (U) 

The minimization of (0 using (TBI) , as done for the [2/2] ap- 
proximant, yields the numerically optimal parameters, 



a = 0.2598, b = 0.5446, c = 1 .9932, a 4 = 0.0538, 
b A = 0.1226, w = 0.1839, A = 3.9439, 



(15) 



The solution (0 with A(x) specified by ( TBI and ( TT3T > will be re- 
ferred as the Hyb(4/4) in this paper and is also plotted in Fig.l. 
The Fig. 2 shows the profiles of F(x) at large x for Hyb(2/2) as 
well as Hyb(4/4), and the numeric solution. The asymptotic ex- 
pansion of the solution Hyb(4/4) shows that for small x <K 1 the 
profile becomes while for x — > oo 

F w (x) =n- 2.0243x + 0.4654x 3 - 0.4270x 5 + • • • . 



while for x — » 

F w (x) 



2.220 



0.9150 9.5908 65.358 



x L 

264.78 



7T 

821.88 



Here, a better value F w (50) = 8.877 x 10~ 4 is obtained for the 
latter asymptotic profile in contrast with the solution Hyb(2/2). 
The computed Skyrmion energies (0, measured in the unit of 
2f„/e, are listed in Table I, including the corresponding results 
obtained by the numeric solution and obtained in the relevant 
references as indicated. 
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In solving (01 numerically, we employ the nonlinear shoot 
algorithm for the boundary values at x — 0.001 and x = 1000 
based on the asymptotic formulas ( ITOb and ( fT2T > of the chiral 
angle F(x). 

3. The static properties of nucleons at low energy 

The static properties of nucleons can be extracted by semi- 
classically quantizing the spinning modes of Skyrme La- 
grangian using the collective variables yy. Here, we will use 
the solution Hyb(2/2) and Hyb(4/4) to compute the static prop- 
erties of nucleons and nucleon-isobar(A) within the framework 
of the bosonic quantization of Skyrme model. 

Following Adkin et al.j^], one can choose S t/(2)-variable 
A(t) as the collective variables, and substitute U — 
A(t)Uo(x)A(t) Jf into (|TJ. In the adiabatic limit, one has 



S SK = S S 



dtTr\ 



r dA <9A + 
~dt~dT 



(16) 



with Sq K the action for the static hedgehog configuration, Iq 

n/(.3e 3 f„), and 



Jr»oo 
x : 
o 



dxsin 2 F[l + F 2 X + sin 2 F/x 2 ]. 



(17) 



which is independent of f„ and e. The Hamiltonian associ- 
ated to (fTol l. when quantized via the quantization procedure 
in terms of collective coordinates, yields an eigenvalue (H) = 
M + J(J + l)/(27oA), with M - E SK being the soli ton energy 
of the Skyrmion. This yields the masses of the nucleon and 
A-isobar 

M /v=M+-^-,Ma = M+-^-. (18) 
8/oA 8/ ( )A 

By adjusting f n and e to fit the nucleon and delta masses in ( fT8l . 
one can fix the model parameters f„ and e using the calculated 
M and A through © and ( TPTT i. 

The isoscalor root mean square(r.m.s) radius and isoscalor 
magnetic r.m.s radius are given by 



efn{r 2 ))L 2 ={— x 2 sin 2 FF X 



1/2 



efir 2 ) 1 ' 2 = 



C x A sin 2 FF x dx 1 
C° x 2 sin 2 FF x dx j 



1/2 



(19) 



respectively. Combining with the masses of nucleon and the 
delta, one can evaluate the magnetic moments for proton and 
neutron via the following formula 



flp,n 



<r 2 >/= 



M N (M A - M N ) : 



2{M A -M N Y 
(20) 

where plus and minus correspond to proton and neutron, respec- 
tively. The calculated results for these quantities using two so- 
lution schemes (Hyb.(2/2) and Hyb.(4/4)) are shown explicitly 
in Table II, compared to the experimental values as well as that 
computed by the numeric solution for F(x). The corresponding 
results from other predictions are also shown in this table. Here 



in Table II, we use the experimental values = 938.9MeV, 
Ma = 1232MeV for fixing e and f n through ( fTSl . in contrast 
with me input Ma, = 938MeV,M A = 1232MeV usedby Ref.Q 
and Ref.lfll. 

To checkthe solution further, we also list in the table, the 
axial coupling constant and the 7rAW-coupling, which are given 
by 

n „ M N 

lJ,gnNN 



gA 3<j2 ^ 

respectively. Here, the numeric factor G is 



gA 



(21) 



Jo 



4 I dxxr 
Jo 



F x + 



sin 2F sin 2F 



2 sin 2 F sin 2 F 
r — F x + — r— sm2F 



(22) 



4. Concluding remarks 

We show that the hybrid form of a kink-like solution and 
that given by the instanton method are suited to approximate 
the exact solution for the hedgehog Skyrmion, when combining 
with Pade approximation. The resulted analytic solution (0 
has two remarkable features: (1) it is simple in the sense that 
it is globally given in whole region; (2) it well approaches the 
asymptotic behavior of the exact solution. We note that the 
further generalization of ©, made by approximating c in © via 
Pade approximation, does not exhibit remarkable improvement, 
particularly in the asymptotic behavior of the chiral angle F (x) 
at infinity. We expect that our solution can be useful in the 
dynamics study of the Skyrmion evolution and interactions. 
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Table I 



Work 


m 


[13] 


16] 


[14] 


[12] 


Hyb(2/2) 


Hyb(4/4) 


Num. 


(2/,/e) 


36.5 


36.47 


36.484 


36.474 


36.47 


36.4638 


36.4638 


36.4616 


M 


1.233 


1.2317 


1.2322 


1.23186 


1.2317 


1.23152 


1.23146 


1.23145 



Table II 



Quantities 


Ref.[3] 


Ref.[14] 


Hyb.(2/2) 


Hyb.(4/4) 


Num. 


Expt. 


M/(2fJe) 


36.5 


36.5 


36.4638 


36.4638 


36.4616 




2f n (MeV) 


129 


130 


128.730 


129.453 


129.260 


186 


e 


5.45 


5.48 


5.4229 


5.4527 


5.4446 




A 


50.9 


52.2 


50.1467 


51.2830 


50.9782 




(r 2 }^ Q (fm) 


0.59 


0.586 


0.5985 


0.5920 


0.5938 


0.72 




0.92 


0.920 


0.9258 


0.9208 


0.9222 


0.81 


Hp 


1.87 




1.8825 


1.8764 


1.8781 


2.79 


Mn 


-1.31 




-1.3209 


-1.3269 


-1.3253 


-1.91 




1.43 




1.4252 


1.4141 


1.4171 


1.46 


8A 


0.61 




0.6332 


0.5992 


0.6081 


1.23 


gnNN 


8.9 




9.2364 


8.6921 


8.8343 


13.5 
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